Electricity/Solar Update
A look at the latest Jedlovec House eletricity usage and solar production from PPL and Enphase.
Load packages
library(tidyverse)
library(lubridate)
library(hms)
library(readxl)
library(onlineBcp)hourly_ppl_pivot <- read_csv("ppl_15mins.csv")Rows: 234720 Columns: 4── Column specification ────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
chr (1): Read Type
dbl (1): kWh
dttm (1): date
time (1): time
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.(hourly_ppl_net <- hourly_ppl_pivot %>%
filter(`Read Type` == "kWh Net"))NALoad Enphase data
full_hist <- read_csv("enphase_production.csv")New names:Rows: 72672 Columns: 3── Column specification ────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
dbl (2): ...1, Energy Produced (Wh)
dttm (1): Date/Time
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.(hourly_production <- full_hist %>%
rename(datetime = `Date/Time`, energy_produced_Wh = `Energy Produced (Wh)`) %>%
mutate(datetime = as.POSIXct(datetime,format="%m/%d/%Y %H:%M",tz=Sys.timezone())) %>%
mutate(date = date(datetime), time = as.hms(format(datetime, format = "%H:%M:%S")), month = month(datetime, label=TRUE), year = year(datetime), day = day(datetime), yday = yday(datetime), monthday = format(datetime, "%m-%d"), wday = wday(datetime, label=TRUE), equinox_day = (yday + 10) %% 365, equinox_group = floor((equinox_day+15)/30)*30)
)NASpot-check Enphase Should see lifetime production by day and by hour
ggplot(hourly_production, aes(datetime, energy_produced_Wh)) +
geom_point()
ggplot(hourly_production, aes(time, energy_produced_Wh)) +
theme(axis.text.x = element_text(angle = 90)) +
geom_point()
Net + Produced = Consumed
# hourly_ppl_net <- hourly_ppl_net %>% mutate(date = as_date(date))
#hourly_ppl_net %>% arrange(desc(date))
#hourly_production %>% arrange(desc(date))
(hourly_electricity <- hourly_ppl_net %>%
inner_join(hourly_production, by = join_by(date,time)) %>%
mutate(consumed_kWh = kWh + energy_produced_Wh/1000, produced_kWh = energy_produced_Wh/1000) %>%
rename(net_kWh = kWh) %>%
select(datetime, date, time, net_kWh, produced_kWh, consumed_kWh) %>%
arrange(desc(datetime)))NACalculate production for first year of solar panels, second year, etc.
hourly_electricity <- hourly_electricity %>%
mutate(solar_year = as_factor(ceiling(interval(as_date(as.POSIXct('04/15/2022',format="%m/%d/%Y",tz=Sys.timezone())),as_date(date)) / years(1) ))) %>%
mutate(yday = yday(date)) %>%
arrange(datetime)
hourly_electricity %>%
group_by(solar_year) %>%
filter(solar_year == 1 | solar_year == 2) %>%
summarize(net_kWh = sum(net_kWh), produced_kWh = sum(produced_kWh), consumed_kWh = sum(consumed_kWh)) NAInteresting! Produced kWh went down by 5%, but consumed kWh went down by 10% despite the fact that we got an electric car! Let’s explore that further.
yday_order <- c(106:366,1:105)
daily_electricity <- hourly_electricity %>%
filter(solar_year == 1 | solar_year == 2) %>%
group_by(solar_year, yday) %>%
summarize(net_kWh = sum(net_kWh), produced_kWh = sum(produced_kWh), consumed_kWh = sum(consumed_kWh), .groups = "drop") %>%
group_by(neworder = cumsum(yday == 106)) %>%
mutate(plot_date = factor(yday, levels = yday_order, ordered = TRUE)) %>%
#mutate(plot_date = factor(yday, levels = yday_order, labels = parse_date_time(yday, orders = "j"), ordered = TRUE)) %>%
arrange(plot_date)
ggplot(daily_electricity, aes(plot_date, consumed_kWh, group=solar_year, color=solar_year)) +
geom_point() +
geom_smooth(span=0.3) +
#scale_x_continuous(breaks = c(1, 32, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335), labels = month.abb) +
theme(axis.text.x = element_text(angle = 90))
Daily totals
daily_electricity <- hourly_electricity %>%
filter(solar_year == 1 | solar_year == 2) %>%
mutate(solar_year_start_year = as.numeric(solar_year) - 1 + 2021) %>%
mutate(solar_year_start_date = paste('04/15/',solar_year_start_year,sep='')) %>%
mutate(solar_day = interval(as_date(as.POSIXct(paste('04/15/',solar_year_start_year),format="%m/%d/%Y",tz=Sys.timezone())),as_date(date)) / days(1)
) %>%
#filter(solar_day < 367) %>%
group_by(date, solar_day, solar_year) %>%
summarize(net_kWh = sum(net_kWh), produced_kWh = sum(produced_kWh), consumed_kWh = sum(consumed_kWh))`summarise()` has grouped output by 'date', 'solar_day'. You can override using the `.groups` argument.ggplot(daily_electricity, aes(solar_day, consumed_kWh, group = solar_year, color = solar_year)) +
geom_point() +
#scale_x_continuous(breaks = c(1, 32, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335), labels = month.abb) +
theme(axis.text.x = element_text(angle = 90)) 
NA
NALook how our electricity consumption got much less predictable, more uneven, after we purchased an EV and installed a level 2 charger in late January/early February 2023!
Also, a couple of outlier high points in Dec 2022 before we got the EV… hmm.
We probably need to filter out the EV charging and compare that separately.
Average consumption based on time of day
#hourly_electricity %>% summarize(min_date = min(date), max_date = max(date))
electricity_by_time <- hourly_electricity %>%
filter(solar_year == 1 | solar_year == 2) %>%
#filter(date <= as.POSIXct('12/31/2022 00:00',format="%m/%d/%Y %H:%M",tz=Sys.timezone()) ) %>%
group_by(time, solar_year) %>%
summarize(produced_kWh = mean(produced_kWh), consumed_kWh = mean(consumed_kWh), net_kWh = mean(net_kWh)) %>%
arrange(time)`summarise()` has grouped output by 'time'. You can override using the `.groups` argument.electricity_by_time
ggplot(electricity_by_time, aes(x=time, y=consumed_kWh, group=solar_year, color=solar_year)) +
geom_point() 
#labs(colour="",x="Time of Day",y="Electricity Consumption (kWh)")+
#scale_color_manual(values = c("red","black","green")) +
#ggtitle("Home Electricity in 15-Minute Intervals (Since April 15, 2022)")
Wow, look how much less was used during the day in year 2!!! Like a 50% reduction!
Is this due to milder weather?
Let’s look at it by month:
#electricity_by_month <-
hourly_electricity %>%
mutate(month = month(date)) %>%
group_by(solar_year, month) %>%
summarize(consumed_kWh = mean(consumed_kWh)) %>%
arrange(month, solar_year) %>%
pivot_wider(names_from = month, values_from = consumed_kWh) %>%
arrange(solar_year)`summarise()` has grouped output by 'solar_year'. You can override using the `.groups` argument.# hourly_electricity %>%
# group_by(solar_day) %>%
# summarize(consumed_kWh = mean(consumed_kWh)) %>%
# arrange(month, solar_year) %>%
# pivot_wider(names_from = month, values_from = consumed_kWh) %>%
# arrange(solar_year)
#Look at limited window during day (no EV charging)
# hourly_electricity %>%
# filter(time == '12:00:00') %>%
# #filter(time > as.POSIXct('04:00:00',format="%H:%M:$s",tz=Sys.timezone()) ) %>% #& time < as.POSIXct('18:00',format="%H:%M",tz=Sys.timezone())) %>%
# mutate(month = month(date)) %>%
# group_by(solar_year, month) %>%
# summarize(consumed_kWh = mean(consumed_kWh)) %>%
# arrange(month, solar_year) %>%
# pivot_wider(names_from = month, values_from = consumed_kWh) %>%
# arrange(solar_year)
Look at specific months
electricity_by_month_time <- hourly_electricity %>%
mutate(month = as_factor(month(date))) %>%
filter(solar_year == 1 | solar_year == 2) %>%
#Sample month
filter(month == 6) %>%
group_by(solar_year, month, time) %>%
summarize(produced_kWh = mean(produced_kWh), consumed_kWh = mean(consumed_kWh), net_kWh = mean(net_kWh)) %>%
arrange(solar_year, month, time)`summarise()` has grouped output by 'solar_year', 'month'. You can override using the `.groups` argument.ggplot(electricity_by_month_time, aes(time, consumed_kWh, group=solar_year, color=solar_year)) +
#facet_grid(rows = vars(month)) +
geom_point(size=0.5) +
ggtitle("Home Electricity in 15-Minute Intervals") +
ylab("Electricity Consumption (kWh)")
Example day with EV charging
sample_day <- hourly_electricity %>%
filter(solar_year == 1 | solar_year == 2) %>%
filter(yday == 161) %>%
group_by(date, solar_year, yday, time) %>%
summarize(produced_kWh = mean(produced_kWh), consumed_kWh = mean(consumed_kWh), net_kWh = mean(net_kWh)) %>%
arrange(date, solar_year, yday, time)`summarise()` has grouped output by 'date', 'solar_year', 'yday'. You can override using the `.groups` argument.ggplot(sample_day, aes(time, consumed_kWh, group=solar_year, color=solar_year)) +
#facet_grid(rows = vars(month)) +
geom_point(size=0.5) +
ggtitle("Home Electricity in 15-Minute Intervals") +
ylab("Electricity Consumption (kWh)")
Ok, let’s detect and remove EV charging sessions.
Test on electricity consumption data
# library(onlineBcp)
#Select a sample to test
sample <- hourly_electricity %>%
mutate(solar_day = interval(as_date(as.POSIXct('04/15/2022',format="%m/%d/%Y",tz=Sys.timezone())),as_date(date)) / days(1) ) %>%
filter(solar_year == 2) %>%
filter(yday < 38 & yday > 35)
ggplot(sample, aes(datetime, consumed_kWh, group=solar_year, color=solar_year)) +
#facet_grid(rows = vars(month)) +
geom_point(size=0.5) +
ggtitle("Home Electricity in 15-Minute Intervals") +
ylab("Electricity Consumption (kWh)")
bcp model comparisons: https://lindeloev.github.io/mcp/articles/packages.html
min_time <- pull(sample %>% summarize(min_time = min(datetime)))
sample <- sample %>% mutate(time_x = interval(min_time,datetime)/minutes(15)+1)
#Filter data to feed to model
x <- sample %>%
ungroup() %>%
select(consumed_kWh)
x <- x$consumed_kWh
#Run online Bayesian Change Point model
#"Online" means in progress, not retroactive
bcp <- online_cp(x)
plot(summary(bcp))[1] "Change points"[1] "Segments"
begin end mean SD LL of CI UL of CI
[1,] 1 22 0.3472727 0.1636542 0.2898818 0.4046636
[2,] 23 71 0.3963061 0.6793674 0.2366690 0.5559433
[3,] 72 86 2.6233333 0.2818983 2.5036113 2.7430554
[4,] 87 118 0.3384375 0.1960597 0.2814289 0.3954461
[5,] 119 131 1.6671538 0.6909474 1.3519434 1.9823643
[6,] 132 192 0.4151148 0.2525212 0.3619333 0.4682962
Not too bad! It picks up the EV charging sessions, but it also picks up some other, presumably HVAC-related consumption changes. It looks like I should be able to weed these out easily.
#Select a sample to test
sample <- hourly_electricity %>%
#filter(solar_year == 2) %>%
mutate(solar_day = interval(as_date(as.POSIXct('04/15/2022',format="%m/%d/%Y",tz=Sys.timezone())),as_date(date)) / days(1) )
ggplot(sample, aes(datetime, consumed_kWh, group=solar_year, color=solar_year)) +
#facet_grid(rows = vars(month)) +
geom_point(size=0.5) +
ggtitle("Home Electricity in 15-Minute Intervals") +
ylab("Electricity Consumption (kWh)")
You can pick out our family vacation in August, and you can also see a period of high electricity usage in February, presumably related to cold weather. - Verify with weather data?
#sample %>% arrange(datetime)
min_time <- pull(sample %>% summarize(min_time = min(datetime)))
sample <- sample %>%
mutate(time_x = interval(min_time,datetime)/minutes(15)+1) %>%
arrange(datetime)
#Filter data to feed to model
x <- sample %>%
ungroup() %>%
select(consumed_kWh)
x <- x$consumed_kWh
#Run online Bayesian Change Point model
#"Online" means in progress, not retroactive
bcp <- online_cp(x)
summary(bcp)[1] "Change points"[1] "Segments"
begin end mean SD LL of CI UL of CI
[1,] 1 68 0.11851471 0.05251481 0.10803969 0.12898972
[2,] 69 89 0.70085714 0.55165300 0.50284875 0.89886553
[3,] 90 233 0.14222917 0.09243684 0.12955874 0.15489959
[4,] 234 240 0.91714286 0.55476224 0.57224926 1.26203645
[5,] 241 408 0.37915476 0.27470555 0.34429375 0.41401578
[6,] 409 410 1.73000000 0.66468037 0.95691880 2.50308120
[7,] 411 524 0.45407895 0.32414913 0.40414226 0.50401563
[8,] 525 526 1.84250000 0.59609102 1.14919420 2.53580580
[9,] 527 599 0.17747945 0.06771880 0.16444253 0.19051638
[10,] 600 602 1.19966667 0.60550007 0.62464948 1.77468385
[11,] 603 696 0.22370213 0.15978499 0.19659401 0.25081024
[12,] 697 724 0.26989286 0.26216920 0.18839805 0.35138766
[13,] 725 792 0.21539706 0.23737175 0.16804904 0.26274508
[14,] 793 843 0.18211765 0.20866132 0.13405761 0.23017769
[15,] 844 888 0.14957778 0.17634812 0.10633717 0.19281839
[16,] 889 938 0.19734000 0.19177406 0.15273001 0.24194999
[17,] 939 943 0.86680000 0.05580502 0.82574978 0.90785022
[18,] 944 1018 0.17470667 0.10118618 0.15548824 0.19392509
[19,] 1019 1032 0.58428571 0.34496364 0.43263776 0.73593367
[20,] 1033 1142 0.18816364 0.08463429 0.17489039 0.20143689
[21,] 1143 1176 0.18558824 0.15582024 0.14163289 0.22954358
[22,] 1177 1178 1.59500000 0.60952605 0.88606809 2.30393191
[23,] 1179 1368 0.32744737 0.28211650 0.29378234 0.36111239
[24,] 1369 1375 0.95057143 0.51726327 0.62899080 1.27215205
[25,] 1376 1466 0.21749451 0.25720279 0.17314567 0.26184334
[26,] 1467 1560 0.21709574 0.22015376 0.17974584 0.25444564
[27,] 1561 1605 0.30462222 0.40168567 0.20612878 0.40311567
[28,] 1606 1714 0.17048624 0.14168952 0.14816328 0.19280920
[29,] 1715 1955 0.14351452 0.11244080 0.13160093 0.15542811
[30,] 1956 1959 0.84625000 0.16573950 0.70994139 0.98255861
[31,] 1960 2067 0.18087037 0.09166977 0.16636123 0.19537951
[32,] 2068 2144 0.19009091 0.11276725 0.16895286 0.21122896
[33,] 2145 2147 1.92800000 0.91683314 1.05732328 2.79867672
[34,] 2148 2232 0.44197647 0.43827136 0.36378473 0.52016821
[35,] 2233 2252 0.79275000 0.43621348 0.63231049 0.95318951
[36,] 2253 2328 0.15798684 0.11320336 0.13662790 0.17934579
[37,] 2329 2335 0.90457143 0.42456050 0.64062376 1.16851909
[38,] 2336 2472 0.20097810 0.18653806 0.17476403 0.22719217
[39,] 2473 2659 0.17288235 0.10114196 0.16071663 0.18504807
[40,] 2660 2670 0.43463636 0.27245266 0.29951564 0.56975709
[41,] 2671 2672 0.33850000 0.18031223 0.12878116 0.54821884
[42,] 2673 2676 0.37275000 0.11853375 0.27526467 0.47023533
[43,] 2677 3132 0.27140132 0.12122497 0.26206368 0.28073895
[44,] 3133 3414 0.24971277 0.10382509 0.23954314 0.25988239
[45,] 3415 3485 0.26271831 0.08814455 0.24551177 0.27992485
[46,] 3486 3548 0.78198413 0.43995383 0.69081158 0.87315667
[47,] 3549 3634 0.51816279 0.27488166 0.46940727 0.56691831
[48,] 3635 3718 0.33130952 0.13270603 0.30749301 0.35512604
[49,] 3719 3721 0.17733333 0.04707795 0.13262544 0.22204122
[50,] 3722 4070 0.15686819 0.10889972 0.14727990 0.16645649
[51,] 4071 4077 0.79457143 0.04656485 0.76562223 0.82352062
[52,] 4078 4292 0.15275349 0.10925063 0.14049796 0.16500902
[53,] 4293 4445 0.38272549 0.31320973 0.34107532 0.42437566
[54,] 4446 4544 0.58247475 0.29179862 0.53423635 0.63071315
[55,] 4545 4548 1.27175000 0.09571268 1.19303333 1.35046667
[56,] 4549 4863 0.39110159 0.18810262 0.37366880 0.40853437
[57,] 4864 4963 0.47486000 0.26307283 0.43158837 0.51813163
[58,] 4964 4980 1.02270588 0.49126759 0.82672175 1.21869001
[59,] 4981 5146 0.35089759 0.13298570 0.33391992 0.36787526
[60,] 5147 5636 0.41125306 0.23248036 0.39397816 0.42852796
[61,] 5637 5718 0.44039024 0.22202345 0.40006108 0.48071940
[62,] 5719 6011 0.49818430 0.28837213 0.47047365 0.52589495
[63,] 6012 6017 1.02733333 0.18042912 0.90617361 1.14849305
[64,] 6018 6230 0.42732864 0.23357544 0.40100387 0.45365341
[65,] 6231 6235 1.04000000 0.11135529 0.95808695 1.12191305
[66,] 6236 6395 0.30019375 0.11013665 0.28587190 0.31451560
[67,] 6396 6558 0.40439877 0.22889281 0.37490939 0.43388815
[68,] 6559 6565 0.82328571 0.09350528 0.76515383 0.88141760
[69,] 6566 6679 0.27366667 0.05711093 0.26486846 0.28246487
[70,] 6680 6683 0.85625000 0.06570832 0.80220972 0.91029028
[71,] 6684 6803 0.28623333 0.12169632 0.26796015 0.30450651
[72,] 6804 7146 0.52422449 0.29807168 0.49775163 0.55069735
[73,] 7147 7350 0.31261275 0.26878429 0.28165880 0.34356669
[74,] 7351 7671 0.42755452 0.29384902 0.40057719 0.45453185
[75,] 7672 8109 0.40544749 0.32093285 0.38022405 0.43067093
[76,] 8110 8312 0.38780788 0.32644096 0.35012155 0.42549422
[77,] 8313 8516 0.40003922 0.33977798 0.36090944 0.43916899
[78,] 8517 8893 0.37877454 0.28384803 0.35472856 0.40282051
[79,] 8894 9536 0.44813997 0.32546790 0.42702794 0.46925200
[80,] 9537 9765 0.65152402 0.32387550 0.61632037 0.68672766
[81,] 9766 10068 0.49356106 0.25400972 0.46955857 0.51756354
[82,] 10069 10134 0.55030303 0.34243544 0.48097099 0.61963507
[83,] 10135 10140 1.35683333 0.22070017 1.20863125 1.50503542
[84,] 10141 10236 0.36204167 0.16197329 0.33485005 0.38923328
[85,] 10237 10332 0.43502083 0.26903739 0.38985560 0.48018607
[86,] 10333 10437 0.35643810 0.17702993 0.32802103 0.38485517
[87,] 10438 10517 0.37961250 0.18511826 0.34556921 0.41365579
[88,] 10518 10591 0.55533784 0.37381261 0.48386100 0.62681468
[89,] 10592 10928 0.59558457 0.32778050 0.56621515 0.62495399
[90,] 10929 10937 1.56755556 0.25295015 1.42886690 1.70624421
[91,] 10938 11783 0.52851655 0.30350088 0.51135320 0.54567990
[92,] 11784 11969 0.49110753 0.26559050 0.45907558 0.52313947
[93,] 11970 12056 0.39251724 0.17475344 0.36170003 0.42333445
[94,] 12057 12131 0.45756000 0.25458680 0.40920599 0.50591401
[95,] 12132 12726 0.51640672 0.26338135 0.49864628 0.53416717
[96,] 12727 12732 1.29900000 0.20395686 1.16204120 1.43595880
[97,] 12733 13382 0.51832462 0.27554900 0.50054717 0.53610206
[98,] 13383 13603 0.43812670 0.24908223 0.41056704 0.46568636
[99,] 13604 13870 0.43283521 0.20554365 0.41214448 0.45352593
[100,] 13871 14255 0.40451948 0.21354805 0.38661785 0.42242111
[101,] 14256 14551 0.29983784 0.21176357 0.27959214 0.32008353
[102,] 14552 14650 0.32707071 0.28589020 0.27980905 0.37433236
[103,] 14651 14656 0.99883333 0.15830656 0.89252910 1.10513756
[104,] 14657 14964 0.31847727 0.12287769 0.30696065 0.32999389
[105,] 14965 15052 0.38471591 0.16570303 0.35566121 0.41377061
[106,] 15053 15059 1.37171429 0.11153731 1.30237196 1.44105661
[107,] 15060 15128 0.41749275 0.21528795 0.37486206 0.46012345
[108,] 15129 15315 0.39243316 0.20651343 0.36759297 0.41727334
[109,] 15316 15342 0.62296296 0.27143245 0.53704042 0.70888551
[110,] 15343 15376 0.38794118 0.13503713 0.34984854 0.42603381
[111,] 15377 15382 0.31166667 0.08084965 0.25737542 0.36595791
[112,] 15383 15386 0.28250000 0.07889867 0.21761162 0.34738838
[113,] 15387 15389 0.29366667 0.08485478 0.21308376 0.37424957
[114,] 15390 15400 0.27827273 0.04756909 0.25468122 0.30186424
[115,] 15401 15832 0.29526157 0.07287639 0.28949428 0.30102887
[116,] 15833 15838 0.94333333 0.04131182 0.91559209 0.97107458
[117,] 15839 16069 0.31194805 0.08370563 0.30288914 0.32100696
[118,] 16070 16280 0.32904265 0.16637795 0.31020262 0.34788269
[119,] 16281 16307 0.41777778 0.29539365 0.32427026 0.51128530
[120,] 16308 16496 0.30491005 0.06815562 0.29675554 0.31306457
[121,] 16497 16539 0.63139535 0.53210535 0.49792320 0.76486750
[122,] 16540 16602 0.51480952 0.28441444 0.45586973 0.57374932
[123,] 16603 16612 0.36000000 0.01632993 0.35150601 0.36849399
[124,] 16613 16663 0.34460784 0.15144307 0.30972663 0.37948906
[125,] 16664 17012 0.30289112 0.08392813 0.29550150 0.31028074
[126,] 17013 17046 0.60958824 0.29497931 0.52637750 0.69279897
[127,] 17047 17301 0.31431765 0.06981058 0.30712683 0.32150847
[128,] 17302 17309 0.62087500 0.30015779 0.44632017 0.79542983
[129,] 17310 17404 0.29954737 0.07243482 0.28732338 0.31177135
[130,] 17405 17429 0.47876000 0.31626707 0.37471739 0.58280261
[131,] 17430 17650 0.30241176 0.07319617 0.29431299 0.31051054
[132,] 17651 17658 1.01375000 0.02825269 0.99731983 1.03018017
[133,] 17659 17908 0.32274000 0.08432301 0.31396790 0.33151210
[134,] 17909 18014 0.36847170 0.12502687 0.34849712 0.38844628
[135,] 18015 18072 0.45996552 0.22579233 0.41119892 0.50873211
[136,] 18073 18074 1.84500000 0.60104076 1.14593721 2.54406279
[137,] 18075 18173 0.49686869 0.27816552 0.45088403 0.54285334
[138,] 18174 18264 0.32750549 0.08028289 0.31366252 0.34134847
[139,] 18265 18360 0.34258333 0.19949234 0.30909312 0.37607354
[140,] 18361 18433 0.38735616 0.20038306 0.34877931 0.42593302
[141,] 18434 18437 1.30000000 0.06377042 1.24755350 1.35244650
[142,] 18438 18608 0.32561404 0.07849521 0.31574051 0.33548756
[143,] 18609 18615 1.01857143 0.12992672 0.93779648 1.09934638
[144,] 18616 18788 0.30564740 0.07546235 0.29621039 0.31508441
[145,] 18789 18899 0.30182883 0.05804275 0.29276704 0.31089061
[146,] 18900 18905 1.08166667 0.04070217 1.05433481 1.10899853
[147,] 18906 19048 0.34625874 0.09462407 0.33324324 0.35927424
[148,] 19049 19050 1.60950000 0.71771338 0.77473678 2.44426322
[149,] 19051 19261 0.32608531 0.10333138 0.31438444 0.33778618
[150,] 19262 19267 0.92983333 0.05562523 0.89248051 0.96718616
[151,] 19268 19519 0.30967063 0.07631779 0.30176289 0.31757838
[152,] 19520 19624 0.29762857 0.08261440 0.28436720 0.31088994
[153,] 19625 19629 0.96340000 0.11106440 0.88170093 1.04509907
[154,] 19630 19845 0.29996296 0.08392853 0.29056984 0.30935609
[155,] 19846 19852 0.88614286 0.10043311 0.82370397 0.94858174
[156,] 19853 20015 0.28745399 0.07412420 0.27790420 0.29700377
[157,] 20016 20017 1.65850000 0.26657926 1.34844509 1.96855491
[158,] 20018 20092 0.42212000 0.17762259 0.38838390 0.45585610
[159,] 20093 20103 0.87218182 0.47699556 0.63561967 1.10874396
[160,] 20104 20296 0.30757513 0.08093623 0.29799235 0.31715791
[161,] 20297 20340 0.60940909 0.39307880 0.51193694 0.70688124
[162,] 20341 20478 0.36581159 0.20908988 0.33653499 0.39508820
[163,] 20479 20569 0.44713187 0.15490586 0.42042184 0.47384190
[164,] 20570 20578 2.04088889 0.63784823 1.69116657 2.39061121
[165,] 20579 20720 0.57975352 0.26651694 0.54296538 0.61654167
[166,] 20721 20726 1.20166667 0.18071156 1.08031729 1.32301605
[ reached getOption("max.print") -- omitted 790 rows ]plot(summary(bcp))[1] "Change points"[1] "Segments"
begin end mean SD LL of CI UL of CI
[1,] 1 68 0.11851471 0.05251481 0.10803969 0.12898972
[2,] 69 89 0.70085714 0.55165300 0.50284875 0.89886553
[3,] 90 233 0.14222917 0.09243684 0.12955874 0.15489959
[4,] 234 240 0.91714286 0.55476224 0.57224926 1.26203645
[5,] 241 408 0.37915476 0.27470555 0.34429375 0.41401578
[6,] 409 410 1.73000000 0.66468037 0.95691880 2.50308120
[7,] 411 524 0.45407895 0.32414913 0.40414226 0.50401563
[8,] 525 526 1.84250000 0.59609102 1.14919420 2.53580580
[9,] 527 599 0.17747945 0.06771880 0.16444253 0.19051638
[10,] 600 602 1.19966667 0.60550007 0.62464948 1.77468385
[11,] 603 696 0.22370213 0.15978499 0.19659401 0.25081024
[12,] 697 724 0.26989286 0.26216920 0.18839805 0.35138766
[13,] 725 792 0.21539706 0.23737175 0.16804904 0.26274508
[14,] 793 843 0.18211765 0.20866132 0.13405761 0.23017769
[15,] 844 888 0.14957778 0.17634812 0.10633717 0.19281839
[16,] 889 938 0.19734000 0.19177406 0.15273001 0.24194999
[17,] 939 943 0.86680000 0.05580502 0.82574978 0.90785022
[18,] 944 1018 0.17470667 0.10118618 0.15548824 0.19392509
[19,] 1019 1032 0.58428571 0.34496364 0.43263776 0.73593367
[20,] 1033 1142 0.18816364 0.08463429 0.17489039 0.20143689
[21,] 1143 1176 0.18558824 0.15582024 0.14163289 0.22954358
[22,] 1177 1178 1.59500000 0.60952605 0.88606809 2.30393191
[23,] 1179 1368 0.32744737 0.28211650 0.29378234 0.36111239
[24,] 1369 1375 0.95057143 0.51726327 0.62899080 1.27215205
[25,] 1376 1466 0.21749451 0.25720279 0.17314567 0.26184334
[26,] 1467 1560 0.21709574 0.22015376 0.17974584 0.25444564
[27,] 1561 1605 0.30462222 0.40168567 0.20612878 0.40311567
[28,] 1606 1714 0.17048624 0.14168952 0.14816328 0.19280920
[29,] 1715 1955 0.14351452 0.11244080 0.13160093 0.15542811
[30,] 1956 1959 0.84625000 0.16573950 0.70994139 0.98255861
[31,] 1960 2067 0.18087037 0.09166977 0.16636123 0.19537951
[32,] 2068 2144 0.19009091 0.11276725 0.16895286 0.21122896
[33,] 2145 2147 1.92800000 0.91683314 1.05732328 2.79867672
[34,] 2148 2232 0.44197647 0.43827136 0.36378473 0.52016821
[35,] 2233 2252 0.79275000 0.43621348 0.63231049 0.95318951
[36,] 2253 2328 0.15798684 0.11320336 0.13662790 0.17934579
[37,] 2329 2335 0.90457143 0.42456050 0.64062376 1.16851909
[38,] 2336 2472 0.20097810 0.18653806 0.17476403 0.22719217
[39,] 2473 2659 0.17288235 0.10114196 0.16071663 0.18504807
[40,] 2660 2670 0.43463636 0.27245266 0.29951564 0.56975709
[41,] 2671 2672 0.33850000 0.18031223 0.12878116 0.54821884
[42,] 2673 2676 0.37275000 0.11853375 0.27526467 0.47023533
[43,] 2677 3132 0.27140132 0.12122497 0.26206368 0.28073895
[44,] 3133 3414 0.24971277 0.10382509 0.23954314 0.25988239
[45,] 3415 3485 0.26271831 0.08814455 0.24551177 0.27992485
[46,] 3486 3548 0.78198413 0.43995383 0.69081158 0.87315667
[47,] 3549 3634 0.51816279 0.27488166 0.46940727 0.56691831
[48,] 3635 3718 0.33130952 0.13270603 0.30749301 0.35512604
[49,] 3719 3721 0.17733333 0.04707795 0.13262544 0.22204122
[50,] 3722 4070 0.15686819 0.10889972 0.14727990 0.16645649
[51,] 4071 4077 0.79457143 0.04656485 0.76562223 0.82352062
[52,] 4078 4292 0.15275349 0.10925063 0.14049796 0.16500902
[53,] 4293 4445 0.38272549 0.31320973 0.34107532 0.42437566
[54,] 4446 4544 0.58247475 0.29179862 0.53423635 0.63071315
[55,] 4545 4548 1.27175000 0.09571268 1.19303333 1.35046667
[56,] 4549 4863 0.39110159 0.18810262 0.37366880 0.40853437
[57,] 4864 4963 0.47486000 0.26307283 0.43158837 0.51813163
[58,] 4964 4980 1.02270588 0.49126759 0.82672175 1.21869001
[59,] 4981 5146 0.35089759 0.13298570 0.33391992 0.36787526
[60,] 5147 5636 0.41125306 0.23248036 0.39397816 0.42852796
[61,] 5637 5718 0.44039024 0.22202345 0.40006108 0.48071940
[62,] 5719 6011 0.49818430 0.28837213 0.47047365 0.52589495
[63,] 6012 6017 1.02733333 0.18042912 0.90617361 1.14849305
[64,] 6018 6230 0.42732864 0.23357544 0.40100387 0.45365341
[65,] 6231 6235 1.04000000 0.11135529 0.95808695 1.12191305
[66,] 6236 6395 0.30019375 0.11013665 0.28587190 0.31451560
[67,] 6396 6558 0.40439877 0.22889281 0.37490939 0.43388815
[68,] 6559 6565 0.82328571 0.09350528 0.76515383 0.88141760
[69,] 6566 6679 0.27366667 0.05711093 0.26486846 0.28246487
[70,] 6680 6683 0.85625000 0.06570832 0.80220972 0.91029028
[71,] 6684 6803 0.28623333 0.12169632 0.26796015 0.30450651
[72,] 6804 7146 0.52422449 0.29807168 0.49775163 0.55069735
[73,] 7147 7350 0.31261275 0.26878429 0.28165880 0.34356669
[74,] 7351 7671 0.42755452 0.29384902 0.40057719 0.45453185
[75,] 7672 8109 0.40544749 0.32093285 0.38022405 0.43067093
[76,] 8110 8312 0.38780788 0.32644096 0.35012155 0.42549422
[77,] 8313 8516 0.40003922 0.33977798 0.36090944 0.43916899
[78,] 8517 8893 0.37877454 0.28384803 0.35472856 0.40282051
[79,] 8894 9536 0.44813997 0.32546790 0.42702794 0.46925200
[80,] 9537 9765 0.65152402 0.32387550 0.61632037 0.68672766
[81,] 9766 10068 0.49356106 0.25400972 0.46955857 0.51756354
[82,] 10069 10134 0.55030303 0.34243544 0.48097099 0.61963507
[83,] 10135 10140 1.35683333 0.22070017 1.20863125 1.50503542
[84,] 10141 10236 0.36204167 0.16197329 0.33485005 0.38923328
[85,] 10237 10332 0.43502083 0.26903739 0.38985560 0.48018607
[86,] 10333 10437 0.35643810 0.17702993 0.32802103 0.38485517
[87,] 10438 10517 0.37961250 0.18511826 0.34556921 0.41365579
[88,] 10518 10591 0.55533784 0.37381261 0.48386100 0.62681468
[89,] 10592 10928 0.59558457 0.32778050 0.56621515 0.62495399
[90,] 10929 10937 1.56755556 0.25295015 1.42886690 1.70624421
[91,] 10938 11783 0.52851655 0.30350088 0.51135320 0.54567990
[92,] 11784 11969 0.49110753 0.26559050 0.45907558 0.52313947
[93,] 11970 12056 0.39251724 0.17475344 0.36170003 0.42333445
[94,] 12057 12131 0.45756000 0.25458680 0.40920599 0.50591401
[95,] 12132 12726 0.51640672 0.26338135 0.49864628 0.53416717
[96,] 12727 12732 1.29900000 0.20395686 1.16204120 1.43595880
[97,] 12733 13382 0.51832462 0.27554900 0.50054717 0.53610206
[98,] 13383 13603 0.43812670 0.24908223 0.41056704 0.46568636
[99,] 13604 13870 0.43283521 0.20554365 0.41214448 0.45352593
[100,] 13871 14255 0.40451948 0.21354805 0.38661785 0.42242111
[101,] 14256 14551 0.29983784 0.21176357 0.27959214 0.32008353
[102,] 14552 14650 0.32707071 0.28589020 0.27980905 0.37433236
[103,] 14651 14656 0.99883333 0.15830656 0.89252910 1.10513756
[104,] 14657 14964 0.31847727 0.12287769 0.30696065 0.32999389
[105,] 14965 15052 0.38471591 0.16570303 0.35566121 0.41377061
[106,] 15053 15059 1.37171429 0.11153731 1.30237196 1.44105661
[107,] 15060 15128 0.41749275 0.21528795 0.37486206 0.46012345
[108,] 15129 15315 0.39243316 0.20651343 0.36759297 0.41727334
[109,] 15316 15342 0.62296296 0.27143245 0.53704042 0.70888551
[110,] 15343 15376 0.38794118 0.13503713 0.34984854 0.42603381
[111,] 15377 15382 0.31166667 0.08084965 0.25737542 0.36595791
[112,] 15383 15386 0.28250000 0.07889867 0.21761162 0.34738838
[113,] 15387 15389 0.29366667 0.08485478 0.21308376 0.37424957
[114,] 15390 15400 0.27827273 0.04756909 0.25468122 0.30186424
[115,] 15401 15832 0.29526157 0.07287639 0.28949428 0.30102887
[116,] 15833 15838 0.94333333 0.04131182 0.91559209 0.97107458
[117,] 15839 16069 0.31194805 0.08370563 0.30288914 0.32100696
[118,] 16070 16280 0.32904265 0.16637795 0.31020262 0.34788269
[119,] 16281 16307 0.41777778 0.29539365 0.32427026 0.51128530
[120,] 16308 16496 0.30491005 0.06815562 0.29675554 0.31306457
[121,] 16497 16539 0.63139535 0.53210535 0.49792320 0.76486750
[122,] 16540 16602 0.51480952 0.28441444 0.45586973 0.57374932
[123,] 16603 16612 0.36000000 0.01632993 0.35150601 0.36849399
[124,] 16613 16663 0.34460784 0.15144307 0.30972663 0.37948906
[125,] 16664 17012 0.30289112 0.08392813 0.29550150 0.31028074
[126,] 17013 17046 0.60958824 0.29497931 0.52637750 0.69279897
[127,] 17047 17301 0.31431765 0.06981058 0.30712683 0.32150847
[128,] 17302 17309 0.62087500 0.30015779 0.44632017 0.79542983
[129,] 17310 17404 0.29954737 0.07243482 0.28732338 0.31177135
[130,] 17405 17429 0.47876000 0.31626707 0.37471739 0.58280261
[131,] 17430 17650 0.30241176 0.07319617 0.29431299 0.31051054
[132,] 17651 17658 1.01375000 0.02825269 0.99731983 1.03018017
[133,] 17659 17908 0.32274000 0.08432301 0.31396790 0.33151210
[134,] 17909 18014 0.36847170 0.12502687 0.34849712 0.38844628
[135,] 18015 18072 0.45996552 0.22579233 0.41119892 0.50873211
[136,] 18073 18074 1.84500000 0.60104076 1.14593721 2.54406279
[137,] 18075 18173 0.49686869 0.27816552 0.45088403 0.54285334
[138,] 18174 18264 0.32750549 0.08028289 0.31366252 0.34134847
[139,] 18265 18360 0.34258333 0.19949234 0.30909312 0.37607354
[140,] 18361 18433 0.38735616 0.20038306 0.34877931 0.42593302
[141,] 18434 18437 1.30000000 0.06377042 1.24755350 1.35244650
[142,] 18438 18608 0.32561404 0.07849521 0.31574051 0.33548756
[143,] 18609 18615 1.01857143 0.12992672 0.93779648 1.09934638
[144,] 18616 18788 0.30564740 0.07546235 0.29621039 0.31508441
[145,] 18789 18899 0.30182883 0.05804275 0.29276704 0.31089061
[146,] 18900 18905 1.08166667 0.04070217 1.05433481 1.10899853
[147,] 18906 19048 0.34625874 0.09462407 0.33324324 0.35927424
[148,] 19049 19050 1.60950000 0.71771338 0.77473678 2.44426322
[149,] 19051 19261 0.32608531 0.10333138 0.31438444 0.33778618
[150,] 19262 19267 0.92983333 0.05562523 0.89248051 0.96718616
[151,] 19268 19519 0.30967063 0.07631779 0.30176289 0.31757838
[152,] 19520 19624 0.29762857 0.08261440 0.28436720 0.31088994
[153,] 19625 19629 0.96340000 0.11106440 0.88170093 1.04509907
[154,] 19630 19845 0.29996296 0.08392853 0.29056984 0.30935609
[155,] 19846 19852 0.88614286 0.10043311 0.82370397 0.94858174
[156,] 19853 20015 0.28745399 0.07412420 0.27790420 0.29700377
[157,] 20016 20017 1.65850000 0.26657926 1.34844509 1.96855491
[158,] 20018 20092 0.42212000 0.17762259 0.38838390 0.45585610
[159,] 20093 20103 0.87218182 0.47699556 0.63561967 1.10874396
[160,] 20104 20296 0.30757513 0.08093623 0.29799235 0.31715791
[161,] 20297 20340 0.60940909 0.39307880 0.51193694 0.70688124
[162,] 20341 20478 0.36581159 0.20908988 0.33653499 0.39508820
[163,] 20479 20569 0.44713187 0.15490586 0.42042184 0.47384190
[164,] 20570 20578 2.04088889 0.63784823 1.69116657 2.39061121
[165,] 20579 20720 0.57975352 0.26651694 0.54296538 0.61654167
[166,] 20721 20726 1.20166667 0.18071156 1.08031729 1.32301605
[ reached getOption("max.print") -- omitted 790 rows ]
Nice! Looks like it still works over a longer period. (An advantage of an onlineBcp model?)
It looks like any interval with a mean above about 2.0 kWh should qualify as an EV charging session. - Do we need to do anything to avoid mixups with high HVAC usage? - Do we need to put an ideal length of charge on the intervals? Maximum variation during the period? - Use post.prob to filter out any intervals?
result <- summary(bcp, norm.test = TRUE)Warning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov testWarning: ties should not be present for the Kolmogorov-Smirnov test[1] "Change points"[1] "Segments"
begin end mean SD LL of CI UL of CI normality p-value
[1,] 1 68 0.11851471 0.05251481 0.10803969 0.12898972 5.558635e-02
[2,] 69 89 0.70085714 0.55165300 0.50284875 0.89886553 2.356642e-01
[3,] 90 233 0.14222917 0.09243684 0.12955874 0.15489959 5.017888e-02
[4,] 234 240 0.91714286 0.55476224 0.57224926 1.26203645 8.598808e-01
[5,] 241 408 0.37915476 0.27470555 0.34429375 0.41401578 2.983556e-05
[6,] 409 410 1.73000000 0.66468037 0.95691880 2.50308120 9.991595e-01
[7,] 411 524 0.45407895 0.32414913 0.40414226 0.50401563 1.754392e-03
[8,] 525 526 1.84250000 0.59609102 1.14919420 2.53580580 9.991595e-01
[9,] 527 599 0.17747945 0.06771880 0.16444253 0.19051638 2.609959e-02
[10,] 600 602 1.19966667 0.60550007 0.62464948 1.77468385 9.460188e-01
[11,] 603 696 0.22370213 0.15978499 0.19659401 0.25081024 2.906650e-04
[12,] 697 724 0.26989286 0.26216920 0.18839805 0.35138766 1.914599e-02
[13,] 725 792 0.21539706 0.23737175 0.16804904 0.26274508 1.569400e-04
[14,] 793 843 0.18211765 0.20866132 0.13405761 0.23017769 3.147528e-03
[15,] 844 888 0.14957778 0.17634812 0.10633717 0.19281839 4.439415e-04
[16,] 889 938 0.19734000 0.19177406 0.15273001 0.24194999 5.771098e-05
[17,] 939 943 0.86680000 0.05580502 0.82574978 0.90785022 5.927433e-01
[18,] 944 1018 0.17470667 0.10118618 0.15548824 0.19392509 1.033167e-01
[19,] 1019 1032 0.58428571 0.34496364 0.43263776 0.73593367 2.792816e-01
[20,] 1033 1142 0.18816364 0.08463429 0.17489039 0.20143689 4.438916e-03
[21,] 1143 1176 0.18558824 0.15582024 0.14163289 0.22954358 1.927116e-02
[22,] 1177 1178 1.59500000 0.60952605 0.88606809 2.30393191 9.991595e-01
[23,] 1179 1368 0.32744737 0.28211650 0.29378234 0.36111239 1.199396e-07
[24,] 1369 1375 0.95057143 0.51726327 0.62899080 1.27215205 3.258514e-01
[25,] 1376 1466 0.21749451 0.25720279 0.17314567 0.26184334 3.567013e-11
[26,] 1467 1560 0.21709574 0.22015376 0.17974584 0.25444564 6.201917e-06
[27,] 1561 1605 0.30462222 0.40168567 0.20612878 0.40311567 1.684396e-03
[28,] 1606 1714 0.17048624 0.14168952 0.14816328 0.19280920 6.067085e-04
[29,] 1715 1955 0.14351452 0.11244080 0.13160093 0.15542811 1.000881e-06
[30,] 1956 1959 0.84625000 0.16573950 0.70994139 0.98255861 7.932280e-01
[31,] 1960 2067 0.18087037 0.09166977 0.16636123 0.19537951 1.412459e-01
[32,] 2068 2144 0.19009091 0.11276725 0.16895286 0.21122896 1.848765e-01
[33,] 2145 2147 1.92800000 0.91683314 1.05732328 2.79867672 7.289916e-01
[34,] 2148 2232 0.44197647 0.43827136 0.36378473 0.52016821 1.045416e-03
[35,] 2233 2252 0.79275000 0.43621348 0.63231049 0.95318951 9.509594e-01
[36,] 2253 2328 0.15798684 0.11320336 0.13662790 0.17934579 3.341529e-03
[37,] 2329 2335 0.90457143 0.42456050 0.64062376 1.16851909 7.828116e-01
[38,] 2336 2472 0.20097810 0.18653806 0.17476403 0.22719217 3.578956e-08
[39,] 2473 2659 0.17288235 0.10114196 0.16071663 0.18504807 4.748881e-04
[40,] 2660 2670 0.43463636 0.27245266 0.29951564 0.56975709 9.083194e-01
[41,] 2671 2672 0.33850000 0.18031223 0.12878116 0.54821884 9.991595e-01
[42,] 2673 2676 0.37275000 0.11853375 0.27526467 0.47023533 9.750718e-01
[43,] 2677 3132 0.27140132 0.12122497 0.26206368 0.28073895 5.356035e-09
[44,] 3133 3414 0.24971277 0.10382509 0.23954314 0.25988239 8.637535e-14
[45,] 3415 3485 0.26271831 0.08814455 0.24551177 0.27992485 9.550424e-07
[46,] 3486 3548 0.78198413 0.43995383 0.69081158 0.87315667 1.218112e-01
[47,] 3549 3634 0.51816279 0.27488166 0.46940727 0.56691831 2.263517e-02
[48,] 3635 3718 0.33130952 0.13270603 0.30749301 0.35512604 1.468338e-02
[49,] 3719 3721 0.17733333 0.04707795 0.13262544 0.22204122 8.058098e-01
[50,] 3722 4070 0.15686819 0.10889972 0.14727990 0.16645649 1.034916e-05
[51,] 4071 4077 0.79457143 0.04656485 0.76562223 0.82352062 9.627227e-01
[52,] 4078 4292 0.15275349 0.10925063 0.14049796 0.16500902 1.564990e-02
[53,] 4293 4445 0.38272549 0.31320973 0.34107532 0.42437566 8.694972e-03
[54,] 4446 4544 0.58247475 0.29179862 0.53423635 0.63071315 1.060064e-02
[55,] 4545 4548 1.27175000 0.09571268 1.19303333 1.35046667 8.235734e-01
[56,] 4549 4863 0.39110159 0.18810262 0.37366880 0.40853437 1.393781e-05
[57,] 4864 4963 0.47486000 0.26307283 0.43158837 0.51813163 6.418883e-04
[58,] 4964 4980 1.02270588 0.49126759 0.82672175 1.21869001 9.594757e-01
[59,] 4981 5146 0.35089759 0.13298570 0.33391992 0.36787526 4.708571e-05
[60,] 5147 5636 0.41125306 0.23248036 0.39397816 0.42852796 1.728617e-13
[61,] 5637 5718 0.44039024 0.22202345 0.40006108 0.48071940 1.462427e-02
[62,] 5719 6011 0.49818430 0.28837213 0.47047365 0.52589495 4.459741e-04
[63,] 6012 6017 1.02733333 0.18042912 0.90617361 1.14849305 6.175601e-01
[64,] 6018 6230 0.42732864 0.23357544 0.40100387 0.45365341 6.726460e-03
[65,] 6231 6235 1.04000000 0.11135529 0.95808695 1.12191305 9.340271e-01
[66,] 6236 6395 0.30019375 0.11013665 0.28587190 0.31451560 7.056873e-06
[67,] 6396 6558 0.40439877 0.22889281 0.37490939 0.43388815 2.036173e-03
[68,] 6559 6565 0.82328571 0.09350528 0.76515383 0.88141760 3.660591e-01
[69,] 6566 6679 0.27366667 0.05711093 0.26486846 0.28246487 6.685227e-02
[70,] 6680 6683 0.85625000 0.06570832 0.80220972 0.91029028 7.615567e-01
[71,] 6684 6803 0.28623333 0.12169632 0.26796015 0.30450651 1.462045e-03
[72,] 6804 7146 0.52422449 0.29807168 0.49775163 0.55069735 1.414514e-05
[73,] 7147 7350 0.31261275 0.26878429 0.28165880 0.34356669 1.471900e-05
[74,] 7351 7671 0.42755452 0.29384902 0.40057719 0.45453185 4.247697e-02
[75,] 7672 8109 0.40544749 0.32093285 0.38022405 0.43067093 1.267526e-04
[76,] 8110 8312 0.38780788 0.32644096 0.35012155 0.42549422 4.003263e-03
[77,] 8313 8516 0.40003922 0.33977798 0.36090944 0.43916899 3.147864e-03
[78,] 8517 8893 0.37877454 0.28384803 0.35472856 0.40282051 9.653318e-04
[79,] 8894 9536 0.44813997 0.32546790 0.42702794 0.46925200 5.326316e-05
[80,] 9537 9765 0.65152402 0.32387550 0.61632037 0.68672766 7.852331e-03
[81,] 9766 10068 0.49356106 0.25400972 0.46955857 0.51756354 1.900428e-02
[82,] 10069 10134 0.55030303 0.34243544 0.48097099 0.61963507 6.081360e-02
[83,] 10135 10140 1.35683333 0.22070017 1.20863125 1.50503542 4.721032e-01
[84,] 10141 10236 0.36204167 0.16197329 0.33485005 0.38923328 4.720172e-03
[85,] 10237 10332 0.43502083 0.26903739 0.38985560 0.48018607 3.120705e-02
[86,] 10333 10437 0.35643810 0.17702993 0.32802103 0.38485517 1.430184e-02
[87,] 10438 10517 0.37961250 0.18511826 0.34556921 0.41365579 7.363339e-02
[88,] 10518 10591 0.55533784 0.37381261 0.48386100 0.62681468 4.075243e-02
[89,] 10592 10928 0.59558457 0.32778050 0.56621515 0.62495399 6.183850e-03
[90,] 10929 10937 1.56755556 0.25295015 1.42886690 1.70624421 8.302274e-01
[91,] 10938 11783 0.52851655 0.30350088 0.51135320 0.54567990 5.097511e-11
[92,] 11784 11969 0.49110753 0.26559050 0.45907558 0.52313947 9.979034e-03
[93,] 11970 12056 0.39251724 0.17475344 0.36170003 0.42333445 7.095404e-02
[94,] 12057 12131 0.45756000 0.25458680 0.40920599 0.50591401 1.908937e-02
[95,] 12132 12726 0.51640672 0.26338135 0.49864628 0.53416717 2.628987e-09
[96,] 12727 12732 1.29900000 0.20395686 1.16204120 1.43595880 9.826903e-01
[97,] 12733 13382 0.51832462 0.27554900 0.50054717 0.53610206 5.502915e-10
[98,] 13383 13603 0.43812670 0.24908223 0.41056704 0.46568636 3.198705e-07
[99,] 13604 13870 0.43283521 0.20554365 0.41214448 0.45352593 9.398370e-06
[100,] 13871 14255 0.40451948 0.21354805 0.38661785 0.42242111 4.476719e-11
[101,] 14256 14551 0.29983784 0.21176357 0.27959214 0.32008353 2.201160e-10
[102,] 14552 14650 0.32707071 0.28589020 0.27980905 0.37433236 8.394495e-05
[103,] 14651 14656 0.99883333 0.15830656 0.89252910 1.10513756 8.826387e-01
[104,] 14657 14964 0.31847727 0.12287769 0.30696065 0.32999389 2.299932e-02
[105,] 14965 15052 0.38471591 0.16570303 0.35566121 0.41377061 2.742575e-05
[106,] 15053 15059 1.37171429 0.11153731 1.30237196 1.44105661 8.768646e-01
[107,] 15060 15128 0.41749275 0.21528795 0.37486206 0.46012345 5.601728e-03
[108,] 15129 15315 0.39243316 0.20651343 0.36759297 0.41727334 6.875991e-06
[109,] 15316 15342 0.62296296 0.27143245 0.53704042 0.70888551 6.122554e-01
[110,] 15343 15376 0.38794118 0.13503713 0.34984854 0.42603381 5.095012e-01
[111,] 15377 15382 0.31166667 0.08084965 0.25737542 0.36595791 6.752876e-01
[112,] 15383 15386 0.28250000 0.07889867 0.21761162 0.34738838 5.126036e-01
[113,] 15387 15389 0.29366667 0.08485478 0.21308376 0.37424957 8.384185e-01
[114,] 15390 15400 0.27827273 0.04756909 0.25468122 0.30186424 2.677207e-01
[115,] 15401 15832 0.29526157 0.07287639 0.28949428 0.30102887 2.749247e-05
[116,] 15833 15838 0.94333333 0.04131182 0.91559209 0.97107458 9.928218e-01
[117,] 15839 16069 0.31194805 0.08370563 0.30288914 0.32100696 1.579487e-05
[118,] 16070 16280 0.32904265 0.16637795 0.31020262 0.34788269 2.123990e-11
[119,] 16281 16307 0.41777778 0.29539365 0.32427026 0.51128530 1.956821e-03
[120,] 16308 16496 0.30491005 0.06815562 0.29675554 0.31306457 3.143062e-02
[121,] 16497 16539 0.63139535 0.53210535 0.49792320 0.76486750 1.454945e-04
[122,] 16540 16602 0.51480952 0.28441444 0.45586973 0.57374932 1.797704e-03
[123,] 16603 16612 0.36000000 0.01632993 0.35150601 0.36849399 8.186212e-01
[124,] 16613 16663 0.34460784 0.15144307 0.30972663 0.37948906 1.685084e-03
[125,] 16664 17012 0.30289112 0.08392813 0.29550150 0.31028074 1.424762e-06
[126,] 17013 17046 0.60958824 0.29497931 0.52637750 0.69279897 2.099508e-01
[127,] 17047 17301 0.31431765 0.06981058 0.30712683 0.32150847 9.057676e-03
[128,] 17302 17309 0.62087500 0.30015779 0.44632017 0.79542983 8.517411e-01
[129,] 17310 17404 0.29954737 0.07243482 0.28732338 0.31177135 1.722243e-04
[130,] 17405 17429 0.47876000 0.31626707 0.37471739 0.58280261 3.673417e-03
[131,] 17430 17650 0.30241176 0.07319617 0.29431299 0.31051054 1.689783e-05
[132,] 17651 17658 1.01375000 0.02825269 0.99731983 1.03018017 9.393540e-01
[133,] 17659 17908 0.32274000 0.08432301 0.31396790 0.33151210 1.444155e-05
[134,] 17909 18014 0.36847170 0.12502687 0.34849712 0.38844628 3.577015e-02
[135,] 18015 18072 0.45996552 0.22579233 0.41119892 0.50873211 2.914935e-02
[136,] 18073 18074 1.84500000 0.60104076 1.14593721 2.54406279 9.991595e-01
[137,] 18075 18173 0.49686869 0.27816552 0.45088403 0.54285334 3.669993e-04
[138,] 18174 18264 0.32750549 0.08028289 0.31366252 0.34134847 2.043566e-02
[139,] 18265 18360 0.34258333 0.19949234 0.30909312 0.37607354 1.058334e-05
[140,] 18361 18433 0.38735616 0.20038306 0.34877931 0.42593302 2.253355e-04
[141,] 18434 18437 1.30000000 0.06377042 1.24755350 1.35244650 7.318359e-01
[142,] 18438 18608 0.32561404 0.07849521 0.31574051 0.33548756 2.529906e-03
[ reached getOption("max.print") -- omitted 814 rows ]bcp_segments <- as_tibble(result$result$segment) %>%
mutate(length = end - begin) %>%
mutate(lag_mean = lag(mean,1)) %>%
mutate(lead_mean = lead(mean,1)) %>%
mutate(adj = (lag_mean+lead_mean)/2)
bcp_segments %>%
filter(mean > 2.0 & mean < 3.0) NA
sample_adj <- sample %>%
left_join(bcp_segments, by = join_by(time_x>=begin,time_x<=end)) %>%
mutate(consumed_kWh_adj = case_when(
mean > 2.0 & length >= 4 & SD < 1.1 ~ consumed_kWh - (mean - adj),
.default = consumed_kWh
)
) %>%
mutate(consumed_kWh_adj = case_when(
consumed_kWh_adj < 0 ~ 0,
.default = consumed_kWh_adj
)
)
sample_adj %>%
mutate(diff = consumed_kWh - consumed_kWh_adj) %>%
group_by(solar_year) %>%
summarize(diff = sum(diff))
ggplot(sample_adj, aes(datetime, consumed_kWh_adj, group=solar_year, color=solar_year)) +
#facet_grid(rows = vars(month)) +
geom_point(size=0.5) +
ggtitle("Home Electricity in 15-Minute Intervals") +
ylab("Electricity Consumption (kWh)")
solar_year
1 595.8770
2 3504.1387
3 321.9421
vs.Â
solar_year
2 4382
3 607
Underestimates
Ignore the bcp model, just pull data on home charging from Emporia
emporia_daily <- read_csv('/Users/benjedlovec/IoniqCharging/emporia_daily.csv') %>%
mutate(date = date(datetime))Rows: 494 Columns: 4── Column specification ────────────────────────────────────────────────────────────────────────────────────────────
Delimiter: ","
dbl (3): kWh, solar_year, vehicle_year
dttm (1): datetime
ℹ Use `spec()` to retrieve the full column specification for this data.
ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.daily_electricity_adj <- daily_electricity %>%
mutate(date = as_date(date)) %>%
left_join(emporia_daily, by = join_by(date)) %>%
mutate(consumed_kWh_adj = consumed_kWh - replace_na(kWh,0)) %>%
rename(solar_year = solar_year.x, charging_kWh = kWh)
daily_electricity_adj %>%
mutate(diff = consumed_kWh - consumed_kWh_adj) %>%
group_by(solar_year) %>%
summarize(consumed_kWh = sum(consumed_kWh), charging_kWh = sum(replace_na(charging_kWh,0)), consumed_kWh_adj = sum(replace_na(consumed_kWh_adj,0)), diff = sum(replace_na(diff,0)))
ggplot(daily_electricity_adj, aes(solar_day, consumed_kWh_adj, group = solar_year, color = solar_year)) +
geom_point() +
#scale_x_continuous(breaks = c(1, 32, 60, 91, 121, 152, 182, 213, 244, 274, 305, 335), labels = month.abb) +
theme(axis.text.x = element_text(angle = 90)) 
NA
NAmonthly_reduction <- daily_electricity_adj %>%
mutate(year_1 = sum(case_when(solar_year == 1 ~ consumed_kWh_adj)), year_2 = sum(case_when(solar_year == 2 ~ consumed_kWh_adj))) %>%
group_by(month = month(date)) %>%
summarize(year_1 = round(sum(replace_na(year_1,0))), year_2 = round(sum(replace_na(year_2,0)))) %>%
mutate(pct = round((year_1-year_2)/year_1,2))
monthly_reduction
daily_electricity_adj %>%
group_by(month = month(date), solar_year) %>%
summarize(consumed_kWh_adj = sum(consumed_kWh_adj))`summarise()` has grouped output by 'month'. You can override using the `.groups` argument.LS0tCnRpdGxlOiAiRWxlY3RyaWNpdHkvU29sYXIgVXBkYXRlIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgpBIGxvb2sgYXQgdGhlIGxhdGVzdCBKZWRsb3ZlYyBIb3VzZSBlbGV0cmljaXR5IHVzYWdlIGFuZCBzb2xhciBwcm9kdWN0aW9uIGZyb20gUFBMIGFuZCBFbnBoYXNlLgoKTG9hZCBwYWNrYWdlcwpgYGB7cn0KbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkobHVicmlkYXRlKQpsaWJyYXJ5KGhtcykKbGlicmFyeShyZWFkeGwpCmxpYnJhcnkob25saW5lQmNwKQpgYGAKCmBgYHtyfQpob3VybHlfcHBsX3Bpdm90IDwtIHJlYWRfY3N2KCJwcGxfMTVtaW5zLmNzdiIpCgooaG91cmx5X3BwbF9uZXQgPC0gaG91cmx5X3BwbF9waXZvdCAlPiUgCiAgZmlsdGVyKGBSZWFkIFR5cGVgID09ICJrV2ggTmV0IikpCgpgYGAKCkxvYWQgRW5waGFzZSBkYXRhCmBgYHtyfQpmdWxsX2hpc3QgPC0gcmVhZF9jc3YoImVucGhhc2VfcHJvZHVjdGlvbi5jc3YiKQoKKGhvdXJseV9wcm9kdWN0aW9uIDwtIGZ1bGxfaGlzdCAlPiUgCiAgcmVuYW1lKGRhdGV0aW1lID0gYERhdGUvVGltZWAsIGVuZXJneV9wcm9kdWNlZF9XaCA9IGBFbmVyZ3kgUHJvZHVjZWQgKFdoKWApICU+JSAKICBtdXRhdGUoZGF0ZXRpbWUgPSBhcy5QT1NJWGN0KGRhdGV0aW1lLGZvcm1hdD0iJW0vJWQvJVkgJUg6JU0iLHR6PVN5cy50aW1lem9uZSgpKSkgJT4lIAogIG11dGF0ZShkYXRlID0gZGF0ZShkYXRldGltZSksIHRpbWUgPSBhcy5obXMoZm9ybWF0KGRhdGV0aW1lLCBmb3JtYXQgPSAiJUg6JU06JVMiKSksIG1vbnRoID0gbW9udGgoZGF0ZXRpbWUsIGxhYmVsPVRSVUUpLCB5ZWFyID0geWVhcihkYXRldGltZSksIGRheSA9IGRheShkYXRldGltZSksIHlkYXkgPSB5ZGF5KGRhdGV0aW1lKSwgbW9udGhkYXkgPSBmb3JtYXQoZGF0ZXRpbWUsICIlbS0lZCIpLCB3ZGF5ID0gd2RheShkYXRldGltZSwgbGFiZWw9VFJVRSksIGVxdWlub3hfZGF5ID0gKHlkYXkgKyAxMCkgJSUgMzY1LCBlcXVpbm94X2dyb3VwID0gZmxvb3IoKGVxdWlub3hfZGF5KzE1KS8zMCkqMzApCikKCmBgYAoKU3BvdC1jaGVjayBFbnBoYXNlClNob3VsZCBzZWUgbGlmZXRpbWUgcHJvZHVjdGlvbiBieSBkYXkgYW5kIGJ5IGhvdXIKYGBge3J9CmdncGxvdChob3VybHlfcHJvZHVjdGlvbiwgYWVzKGRhdGV0aW1lLCBlbmVyZ3lfcHJvZHVjZWRfV2gpKSArCiAgZ2VvbV9wb2ludCgpCgpnZ3Bsb3QoaG91cmx5X3Byb2R1Y3Rpb24sIGFlcyh0aW1lLCBlbmVyZ3lfcHJvZHVjZWRfV2gpKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCkpICsKICBnZW9tX3BvaW50KCkKCmBgYAoKTmV0ICsgUHJvZHVjZWQgPSBDb25zdW1lZAoKYGBge3J9CiMgaG91cmx5X3BwbF9uZXQgPC0gaG91cmx5X3BwbF9uZXQgJT4lIG11dGF0ZShkYXRlID0gYXNfZGF0ZShkYXRlKSkKCiNob3VybHlfcHBsX25ldCAlPiUgYXJyYW5nZShkZXNjKGRhdGUpKQojaG91cmx5X3Byb2R1Y3Rpb24gJT4lIGFycmFuZ2UoZGVzYyhkYXRlKSkKCihob3VybHlfZWxlY3RyaWNpdHkgPC0gaG91cmx5X3BwbF9uZXQgJT4lIAogICAgaW5uZXJfam9pbihob3VybHlfcHJvZHVjdGlvbiwgYnkgPSBqb2luX2J5KGRhdGUsdGltZSkpICAlPiUgCiAgICBtdXRhdGUoY29uc3VtZWRfa1doID0ga1doICsgZW5lcmd5X3Byb2R1Y2VkX1doLzEwMDAsIHByb2R1Y2VkX2tXaCA9IGVuZXJneV9wcm9kdWNlZF9XaC8xMDAwKSAgJT4lIAogICAgcmVuYW1lKG5ldF9rV2ggPSBrV2gpICU+JSAKICAgIHNlbGVjdChkYXRldGltZSwgZGF0ZSwgdGltZSwgbmV0X2tXaCwgcHJvZHVjZWRfa1doLCBjb25zdW1lZF9rV2gpICU+JQogICAgYXJyYW5nZShkZXNjKGRhdGV0aW1lKSkpCgpgYGAKCkNhbGN1bGF0ZSBwcm9kdWN0aW9uIGZvciBmaXJzdCB5ZWFyIG9mIHNvbGFyIHBhbmVscywgc2Vjb25kIHllYXIsIGV0Yy4KYGBge3J9Cgpob3VybHlfZWxlY3RyaWNpdHkgPC0gaG91cmx5X2VsZWN0cmljaXR5ICU+JSAKICBtdXRhdGUoc29sYXJfeWVhciA9IGFzX2ZhY3RvcihjZWlsaW5nKGludGVydmFsKGFzX2RhdGUoYXMuUE9TSVhjdCgnMDQvMTUvMjAyMicsZm9ybWF0PSIlbS8lZC8lWSIsdHo9U3lzLnRpbWV6b25lKCkpKSxhc19kYXRlKGRhdGUpKSAvIHllYXJzKDEpICkpKSAlPiUgCiAgbXV0YXRlKHlkYXkgPSB5ZGF5KGRhdGUpKSAlPiUgCiAgYXJyYW5nZShkYXRldGltZSkKCmhvdXJseV9lbGVjdHJpY2l0eSAlPiUgCiAgZ3JvdXBfYnkoc29sYXJfeWVhcikgJT4lIAogIGZpbHRlcihzb2xhcl95ZWFyID09IDEgfCBzb2xhcl95ZWFyID09IDIpICU+JSAKICBzdW1tYXJpemUobmV0X2tXaCA9IHN1bShuZXRfa1doKSwgcHJvZHVjZWRfa1doID0gc3VtKHByb2R1Y2VkX2tXaCksIGNvbnN1bWVkX2tXaCA9IHN1bShjb25zdW1lZF9rV2gpKSAKCmBgYAoKSW50ZXJlc3RpbmchIFByb2R1Y2VkIGtXaCB3ZW50IGRvd24gYnkgNSUsIGJ1dCBjb25zdW1lZCBrV2ggd2VudCBkb3duIGJ5IDEwJSBkZXNwaXRlIHRoZSBmYWN0IHRoYXQgd2UgZ290IGFuIGVsZWN0cmljIGNhciEgTGV0J3MgZXhwbG9yZSB0aGF0IGZ1cnRoZXIuIAoKYGBge3J9CnlkYXlfb3JkZXIgPC0gYygxMDY6MzY2LDE6MTA1KQoKZGFpbHlfZWxlY3RyaWNpdHkgPC0gaG91cmx5X2VsZWN0cmljaXR5ICU+JSAKICBmaWx0ZXIoc29sYXJfeWVhciA9PSAxIHwgc29sYXJfeWVhciA9PSAyKSAlPiUgCiAgZ3JvdXBfYnkoc29sYXJfeWVhciwgeWRheSkgJT4lIAogIHN1bW1hcml6ZShuZXRfa1doID0gc3VtKG5ldF9rV2gpLCBwcm9kdWNlZF9rV2ggPSBzdW0ocHJvZHVjZWRfa1doKSwgY29uc3VtZWRfa1doID0gc3VtKGNvbnN1bWVkX2tXaCksIC5ncm91cHMgPSAiZHJvcCIpICU+JSAKICBncm91cF9ieShuZXdvcmRlciA9IGN1bXN1bSh5ZGF5ID09IDEwNikpICU+JSAKICBtdXRhdGUocGxvdF9kYXRlID0gZmFjdG9yKHlkYXksIGxldmVscyA9IHlkYXlfb3JkZXIsIG9yZGVyZWQgPSBUUlVFKSkgJT4lCiAgI211dGF0ZShwbG90X2RhdGUgPSBmYWN0b3IoeWRheSwgbGV2ZWxzID0geWRheV9vcmRlciwgbGFiZWxzID0gcGFyc2VfZGF0ZV90aW1lKHlkYXksIG9yZGVycyA9ICJqIiksIG9yZGVyZWQgPSBUUlVFKSkgJT4lCiAgYXJyYW5nZShwbG90X2RhdGUpCgpnZ3Bsb3QoZGFpbHlfZWxlY3RyaWNpdHksIGFlcyhwbG90X2RhdGUsIGNvbnN1bWVkX2tXaCwgZ3JvdXA9c29sYXJfeWVhciwgY29sb3I9c29sYXJfeWVhcikpICsgCiAgZ2VvbV9wb2ludCgpICsKICBnZW9tX3Ntb290aChzcGFuPTAuMykgKwogICNzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzID0gYygxLCAzMiwgNjAsIDkxLCAxMjEsIDE1MiwgMTgyLCAyMTMsIDI0NCwgMjc0LCAzMDUsIDMzNSksIGxhYmVscyA9IG1vbnRoLmFiYikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTApKSAKCmBgYAoKRGFpbHkgdG90YWxzCmBgYHtyfQoKZGFpbHlfZWxlY3RyaWNpdHkgPC0gaG91cmx5X2VsZWN0cmljaXR5ICU+JSAKICBmaWx0ZXIoc29sYXJfeWVhciA9PSAxIHwgc29sYXJfeWVhciA9PSAyKSAlPiUgCiAgbXV0YXRlKHNvbGFyX3llYXJfc3RhcnRfeWVhciA9IGFzLm51bWVyaWMoc29sYXJfeWVhcikgLSAxICsgMjAyMSkgJT4lIAogIG11dGF0ZShzb2xhcl95ZWFyX3N0YXJ0X2RhdGUgPSBwYXN0ZSgnMDQvMTUvJyxzb2xhcl95ZWFyX3N0YXJ0X3llYXIsc2VwPScnKSkgJT4lIAogIG11dGF0ZShzb2xhcl9kYXkgPSBpbnRlcnZhbChhc19kYXRlKGFzLlBPU0lYY3QocGFzdGUoJzA0LzE1Lycsc29sYXJfeWVhcl9zdGFydF95ZWFyKSxmb3JtYXQ9IiVtLyVkLyVZIix0ej1TeXMudGltZXpvbmUoKSkpLGFzX2RhdGUoZGF0ZSkpIC8gZGF5cygxKSAKICAgICAgICAgKSAlPiUgCiAgI2ZpbHRlcihzb2xhcl9kYXkgPCAzNjcpICU+JSAKICBncm91cF9ieShkYXRlLCBzb2xhcl9kYXksIHNvbGFyX3llYXIpICU+JSAKICBzdW1tYXJpemUobmV0X2tXaCA9IHN1bShuZXRfa1doKSwgcHJvZHVjZWRfa1doID0gc3VtKHByb2R1Y2VkX2tXaCksIGNvbnN1bWVkX2tXaCA9IHN1bShjb25zdW1lZF9rV2gpKQoKZ2dwbG90KGRhaWx5X2VsZWN0cmljaXR5LCBhZXMoc29sYXJfZGF5LCBjb25zdW1lZF9rV2gsIGdyb3VwID0gc29sYXJfeWVhciwgY29sb3IgPSBzb2xhcl95ZWFyKSkgKyAKICBnZW9tX3BvaW50KCkgKwogICNzY2FsZV94X2NvbnRpbnVvdXMoYnJlYWtzID0gYygxLCAzMiwgNjAsIDkxLCAxMjEsIDE1MiwgMTgyLCAyMTMsIDI0NCwgMjc0LCAzMDUsIDMzNSksIGxhYmVscyA9IG1vbnRoLmFiYikgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTApKSAKCiAKYGBgCkxvb2sgaG93IG91ciBlbGVjdHJpY2l0eSBjb25zdW1wdGlvbiBnb3QgbXVjaCBsZXNzIHByZWRpY3RhYmxlLCBtb3JlIHVuZXZlbiwgYWZ0ZXIgd2UgcHVyY2hhc2VkIGFuIEVWIGFuZCBpbnN0YWxsZWQgYSBsZXZlbCAyIGNoYXJnZXIgaW4gbGF0ZSBKYW51YXJ5L2Vhcmx5IEZlYnJ1YXJ5IDIwMjMhCgpBbHNvLCBhIGNvdXBsZSBvZiBvdXRsaWVyIGhpZ2ggcG9pbnRzIGluIERlYyAyMDIyIGJlZm9yZSB3ZSBnb3QgdGhlIEVWLi4uIGhtbS4KCldlIHByb2JhYmx5IG5lZWQgdG8gZmlsdGVyIG91dCB0aGUgRVYgY2hhcmdpbmcgYW5kIGNvbXBhcmUgdGhhdCBzZXBhcmF0ZWx5LiAKCkF2ZXJhZ2UgY29uc3VtcHRpb24gYmFzZWQgb24gdGltZSBvZiBkYXkKYGBge3J9CiNob3VybHlfZWxlY3RyaWNpdHkgJT4lIHN1bW1hcml6ZShtaW5fZGF0ZSA9IG1pbihkYXRlKSwgbWF4X2RhdGUgPSBtYXgoZGF0ZSkpCgplbGVjdHJpY2l0eV9ieV90aW1lIDwtIGhvdXJseV9lbGVjdHJpY2l0eSAlPiUgCiAgZmlsdGVyKHNvbGFyX3llYXIgPT0gMSB8IHNvbGFyX3llYXIgPT0gMikgJT4lIAogICNmaWx0ZXIoZGF0ZSA8PSBhcy5QT1NJWGN0KCcxMi8zMS8yMDIyIDAwOjAwJyxmb3JtYXQ9IiVtLyVkLyVZICVIOiVNIix0ej1TeXMudGltZXpvbmUoKSkgKSAlPiUgCiAgZ3JvdXBfYnkodGltZSwgc29sYXJfeWVhcikgJT4lIAogIHN1bW1hcml6ZShwcm9kdWNlZF9rV2ggPSBtZWFuKHByb2R1Y2VkX2tXaCksIGNvbnN1bWVkX2tXaCA9IG1lYW4oY29uc3VtZWRfa1doKSwgbmV0X2tXaCA9IG1lYW4obmV0X2tXaCkpICU+JSAKICBhcnJhbmdlKHRpbWUpCgplbGVjdHJpY2l0eV9ieV90aW1lCgpnZ3Bsb3QoZWxlY3RyaWNpdHlfYnlfdGltZSwgYWVzKHg9dGltZSwgeT1jb25zdW1lZF9rV2gsIGdyb3VwPXNvbGFyX3llYXIsIGNvbG9yPXNvbGFyX3llYXIpKSArCiAgZ2VvbV9wb2ludCgpIAogICNsYWJzKGNvbG91cj0iIix4PSJUaW1lIG9mIERheSIseT0iRWxlY3RyaWNpdHkgQ29uc3VtcHRpb24gKGtXaCkiKSsKICAjc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGMoInJlZCIsImJsYWNrIiwiZ3JlZW4iKSkgKwogICNnZ3RpdGxlKCJIb21lIEVsZWN0cmljaXR5IGluIDE1LU1pbnV0ZSBJbnRlcnZhbHMgKFNpbmNlIEFwcmlsIDE1LCAyMDIyKSIpCgpgYGAKV293LCBsb29rIGhvdyBtdWNoIGxlc3Mgd2FzIHVzZWQgZHVyaW5nIHRoZSBkYXkgaW4geWVhciAyISEhIExpa2UgYSA1MCUgcmVkdWN0aW9uISAKCklzIHRoaXMgZHVlIHRvIG1pbGRlciB3ZWF0aGVyPyAKCkxldCdzIGxvb2sgYXQgaXQgYnkgbW9udGg6IAoKYGBge3J9CiNlbGVjdHJpY2l0eV9ieV9tb250aCA8LSAKaG91cmx5X2VsZWN0cmljaXR5ICU+JSAKICBtdXRhdGUobW9udGggPSBtb250aChkYXRlKSkgJT4lIAogIGdyb3VwX2J5KHNvbGFyX3llYXIsIG1vbnRoKSAlPiUgCiAgc3VtbWFyaXplKGNvbnN1bWVkX2tXaCA9IG1lYW4oY29uc3VtZWRfa1doKSkgJT4lIAogIGFycmFuZ2UobW9udGgsIHNvbGFyX3llYXIpICU+JSAKICBwaXZvdF93aWRlcihuYW1lc19mcm9tID0gbW9udGgsIHZhbHVlc19mcm9tID0gY29uc3VtZWRfa1doKSAlPiUgCiAgYXJyYW5nZShzb2xhcl95ZWFyKQoKIyBob3VybHlfZWxlY3RyaWNpdHkgJT4lIAojICAgZ3JvdXBfYnkoc29sYXJfZGF5KSAlPiUgCiMgICBzdW1tYXJpemUoY29uc3VtZWRfa1doID0gbWVhbihjb25zdW1lZF9rV2gpKSAlPiUgCiMgICBhcnJhbmdlKG1vbnRoLCBzb2xhcl95ZWFyKSAlPiUgCiMgICBwaXZvdF93aWRlcihuYW1lc19mcm9tID0gbW9udGgsIHZhbHVlc19mcm9tID0gY29uc3VtZWRfa1doKSAlPiUgCiMgICBhcnJhbmdlKHNvbGFyX3llYXIpCgoKCiNMb29rIGF0IGxpbWl0ZWQgd2luZG93IGR1cmluZyBkYXkgKG5vIEVWIGNoYXJnaW5nKQogCiMgaG91cmx5X2VsZWN0cmljaXR5ICU+JSAKIyAgIGZpbHRlcih0aW1lID09ICcxMjowMDowMCcpICU+JSAKIyAgICNmaWx0ZXIodGltZSA+IGFzLlBPU0lYY3QoJzA0OjAwOjAwJyxmb3JtYXQ9IiVIOiVNOiRzIix0ej1TeXMudGltZXpvbmUoKSkgKSAlPiUgICMmIHRpbWUgPCBhcy5QT1NJWGN0KCcxODowMCcsZm9ybWF0PSIlSDolTSIsdHo9U3lzLnRpbWV6b25lKCkpKSAlPiUgCiMgICBtdXRhdGUobW9udGggPSBtb250aChkYXRlKSkgJT4lIAojICAgZ3JvdXBfYnkoc29sYXJfeWVhciwgbW9udGgpICU+JSAKIyAgIHN1bW1hcml6ZShjb25zdW1lZF9rV2ggPSBtZWFuKGNvbnN1bWVkX2tXaCkpICU+JSAKIyAgIGFycmFuZ2UobW9udGgsIHNvbGFyX3llYXIpICU+JSAKIyAgIHBpdm90X3dpZGVyKG5hbWVzX2Zyb20gPSBtb250aCwgdmFsdWVzX2Zyb20gPSBjb25zdW1lZF9rV2gpICU+JSAKIyAgIGFycmFuZ2Uoc29sYXJfeWVhcikKCgpgYGAKCkxvb2sgYXQgc3BlY2lmaWMgbW9udGhzCgpgYGB7cn0KZWxlY3RyaWNpdHlfYnlfbW9udGhfdGltZSA8LSBob3VybHlfZWxlY3RyaWNpdHkgJT4lIAogIG11dGF0ZShtb250aCA9IGFzX2ZhY3Rvcihtb250aChkYXRlKSkpICU+JSAKICBmaWx0ZXIoc29sYXJfeWVhciA9PSAxIHwgc29sYXJfeWVhciA9PSAyKSAlPiUgCiAgI1NhbXBsZSBtb250aAogIGZpbHRlcihtb250aCA9PSA2KSAlPiUgCiAgZ3JvdXBfYnkoc29sYXJfeWVhciwgbW9udGgsIHRpbWUpICU+JSAKICBzdW1tYXJpemUocHJvZHVjZWRfa1doID0gbWVhbihwcm9kdWNlZF9rV2gpLCBjb25zdW1lZF9rV2ggPSBtZWFuKGNvbnN1bWVkX2tXaCksIG5ldF9rV2ggPSBtZWFuKG5ldF9rV2gpKSAlPiUgCiAgYXJyYW5nZShzb2xhcl95ZWFyLCBtb250aCwgdGltZSkKCmdncGxvdChlbGVjdHJpY2l0eV9ieV9tb250aF90aW1lLCBhZXModGltZSwgY29uc3VtZWRfa1doLCBncm91cD1zb2xhcl95ZWFyLCBjb2xvcj1zb2xhcl95ZWFyKSkgKwogICNmYWNldF9ncmlkKHJvd3MgPSB2YXJzKG1vbnRoKSkgKwogIGdlb21fcG9pbnQoc2l6ZT0wLjUpICsKICBnZ3RpdGxlKCJIb21lIEVsZWN0cmljaXR5IGluIDE1LU1pbnV0ZSBJbnRlcnZhbHMiKSArCiAgeWxhYigiRWxlY3RyaWNpdHkgQ29uc3VtcHRpb24gKGtXaCkiKQoKYGBgCgpFeGFtcGxlIGRheSB3aXRoIEVWIGNoYXJnaW5nCgpgYGB7cn0Kc2FtcGxlX2RheSA8LSBob3VybHlfZWxlY3RyaWNpdHkgJT4lIAogIGZpbHRlcihzb2xhcl95ZWFyID09IDEgfCBzb2xhcl95ZWFyID09IDIpICU+JSAKICBmaWx0ZXIoeWRheSA9PSAxNjEpICU+JSAKICBncm91cF9ieShkYXRlLCBzb2xhcl95ZWFyLCB5ZGF5LCB0aW1lKSAlPiUgCiAgc3VtbWFyaXplKHByb2R1Y2VkX2tXaCA9IG1lYW4ocHJvZHVjZWRfa1doKSwgY29uc3VtZWRfa1doID0gbWVhbihjb25zdW1lZF9rV2gpLCBuZXRfa1doID0gbWVhbihuZXRfa1doKSkgJT4lIAogIGFycmFuZ2UoZGF0ZSwgc29sYXJfeWVhciwgeWRheSwgdGltZSkKCmdncGxvdChzYW1wbGVfZGF5LCBhZXModGltZSwgY29uc3VtZWRfa1doLCBncm91cD1zb2xhcl95ZWFyLCBjb2xvcj1zb2xhcl95ZWFyKSkgKwogICNmYWNldF9ncmlkKHJvd3MgPSB2YXJzKG1vbnRoKSkgKwogIGdlb21fcG9pbnQoc2l6ZT0wLjUpICsKICBnZ3RpdGxlKCJIb21lIEVsZWN0cmljaXR5IGluIDE1LU1pbnV0ZSBJbnRlcnZhbHMiKSArCiAgeWxhYigiRWxlY3RyaWNpdHkgQ29uc3VtcHRpb24gKGtXaCkiKQoKYGBgCgpPaywgbGV0J3MgZGV0ZWN0IGFuZCByZW1vdmUgRVYgY2hhcmdpbmcgc2Vzc2lvbnMuIAoKVGVzdCBvbiBlbGVjdHJpY2l0eSBjb25zdW1wdGlvbiBkYXRhCmBgYHtyfQojIGxpYnJhcnkob25saW5lQmNwKQoKI1NlbGVjdCBhIHNhbXBsZSB0byB0ZXN0CnNhbXBsZSA8LSBob3VybHlfZWxlY3RyaWNpdHkgJT4lIAogIG11dGF0ZShzb2xhcl9kYXkgPSBpbnRlcnZhbChhc19kYXRlKGFzLlBPU0lYY3QoJzA0LzE1LzIwMjInLGZvcm1hdD0iJW0vJWQvJVkiLHR6PVN5cy50aW1lem9uZSgpKSksYXNfZGF0ZShkYXRlKSkgLyBkYXlzKDEpICkgJT4lCiAgZmlsdGVyKHNvbGFyX3llYXIgPT0gMikgJT4lIAogIGZpbHRlcih5ZGF5IDwgMzggJiB5ZGF5ID4gMzUpCgpnZ3Bsb3Qoc2FtcGxlLCBhZXMoZGF0ZXRpbWUsIGNvbnN1bWVkX2tXaCwgZ3JvdXA9c29sYXJfeWVhciwgY29sb3I9c29sYXJfeWVhcikpICsKICAjZmFjZXRfZ3JpZChyb3dzID0gdmFycyhtb250aCkpICsKICBnZW9tX3BvaW50KHNpemU9MC41KSArCiAgZ2d0aXRsZSgiSG9tZSBFbGVjdHJpY2l0eSBpbiAxNS1NaW51dGUgSW50ZXJ2YWxzIikgKwogIHlsYWIoIkVsZWN0cmljaXR5IENvbnN1bXB0aW9uIChrV2gpIikKYGBgCgpiY3AgbW9kZWwgY29tcGFyaXNvbnM6IGh0dHBzOi8vbGluZGVsb2V2LmdpdGh1Yi5pby9tY3AvYXJ0aWNsZXMvcGFja2FnZXMuaHRtbAoKCmBgYHtyfQptaW5fdGltZSA8LSBwdWxsKHNhbXBsZSAlPiUgc3VtbWFyaXplKG1pbl90aW1lID0gbWluKGRhdGV0aW1lKSkpCgpzYW1wbGUgPC0gc2FtcGxlICU+JSBtdXRhdGUodGltZV94ID0gaW50ZXJ2YWwobWluX3RpbWUsZGF0ZXRpbWUpL21pbnV0ZXMoMTUpKzEpCgojRmlsdGVyIGRhdGEgdG8gZmVlZCB0byBtb2RlbAp4IDwtIHNhbXBsZSAlPiUgIAogIHVuZ3JvdXAoKSAlPiUgCiAgc2VsZWN0KGNvbnN1bWVkX2tXaCkKCnggPC0geCRjb25zdW1lZF9rV2gKCiNSdW4gb25saW5lIEJheWVzaWFuIENoYW5nZSBQb2ludCBtb2RlbCAKIyJPbmxpbmUiIG1lYW5zIGluIHByb2dyZXNzLCBub3QgcmV0cm9hY3RpdmUKYmNwIDwtIG9ubGluZV9jcCh4KQoKcGxvdChzdW1tYXJ5KGJjcCkpCgpgYGAKTm90IHRvbyBiYWQhIEl0IHBpY2tzIHVwIHRoZSBFViBjaGFyZ2luZyBzZXNzaW9ucywgYnV0IGl0IGFsc28gcGlja3MgdXAgc29tZSBvdGhlciwgcHJlc3VtYWJseSBIVkFDLXJlbGF0ZWQgY29uc3VtcHRpb24gY2hhbmdlcy4gSXQgbG9va3MgbGlrZSBJIHNob3VsZCBiZSBhYmxlIHRvIHdlZWQgdGhlc2Ugb3V0IGVhc2lseS4gCgoKYGBge3J9CgojU2VsZWN0IGEgc2FtcGxlIHRvIHRlc3QKc2FtcGxlIDwtIGhvdXJseV9lbGVjdHJpY2l0eSAlPiUgCiAgI2ZpbHRlcihzb2xhcl95ZWFyID09IDIpICU+JSAKICBtdXRhdGUoc29sYXJfZGF5ID0gaW50ZXJ2YWwoYXNfZGF0ZShhcy5QT1NJWGN0KCcwNC8xNS8yMDIyJyxmb3JtYXQ9IiVtLyVkLyVZIix0ej1TeXMudGltZXpvbmUoKSkpLGFzX2RhdGUoZGF0ZSkpIC8gZGF5cygxKSApIAoKZ2dwbG90KHNhbXBsZSwgYWVzKGRhdGV0aW1lLCBjb25zdW1lZF9rV2gsIGdyb3VwPXNvbGFyX3llYXIsIGNvbG9yPXNvbGFyX3llYXIpKSArCiAgI2ZhY2V0X2dyaWQocm93cyA9IHZhcnMobW9udGgpKSArCiAgZ2VvbV9wb2ludChzaXplPTAuNSkgKwogIGdndGl0bGUoIkhvbWUgRWxlY3RyaWNpdHkgaW4gMTUtTWludXRlIEludGVydmFscyIpICsKICB5bGFiKCJFbGVjdHJpY2l0eSBDb25zdW1wdGlvbiAoa1doKSIpCmBgYApZb3UgY2FuIHBpY2sgb3V0IG91ciBmYW1pbHkgdmFjYXRpb24gaW4gQXVndXN0LCBhbmQgeW91IGNhbiBhbHNvIHNlZSBhIHBlcmlvZCBvZiBoaWdoIGVsZWN0cmljaXR5IHVzYWdlIGluIEZlYnJ1YXJ5LCBwcmVzdW1hYmx5IHJlbGF0ZWQgdG8gY29sZCB3ZWF0aGVyLiAKLSBWZXJpZnkgd2l0aCB3ZWF0aGVyIGRhdGE/IAoKCmBgYHtyfQojc2FtcGxlICU+JSBhcnJhbmdlKGRhdGV0aW1lKQoKbWluX3RpbWUgPC0gcHVsbChzYW1wbGUgJT4lIHN1bW1hcml6ZShtaW5fdGltZSA9IG1pbihkYXRldGltZSkpKQoKc2FtcGxlIDwtIHNhbXBsZSAlPiUgCiAgbXV0YXRlKHRpbWVfeCA9IGludGVydmFsKG1pbl90aW1lLGRhdGV0aW1lKS9taW51dGVzKDE1KSsxKSAlPiUgCiAgYXJyYW5nZShkYXRldGltZSkKCiNGaWx0ZXIgZGF0YSB0byBmZWVkIHRvIG1vZGVsCnggPC0gc2FtcGxlICU+JSAgCiAgdW5ncm91cCgpICU+JSAKICBzZWxlY3QoY29uc3VtZWRfa1doKQoKeCA8LSB4JGNvbnN1bWVkX2tXaAoKI1J1biBvbmxpbmUgQmF5ZXNpYW4gQ2hhbmdlIFBvaW50IG1vZGVsIAojIk9ubGluZSIgbWVhbnMgaW4gcHJvZ3Jlc3MsIG5vdCByZXRyb2FjdGl2ZQpiY3AgPC0gb25saW5lX2NwKHgpCgpzdW1tYXJ5KGJjcCkKCnBsb3Qoc3VtbWFyeShiY3ApKQpgYGAKTmljZSEgTG9va3MgbGlrZSBpdCBzdGlsbCB3b3JrcyBvdmVyIGEgbG9uZ2VyIHBlcmlvZC4gKEFuIGFkdmFudGFnZSBvZiBhbiBvbmxpbmVCY3AgbW9kZWw/KQoKSXQgbG9va3MgbGlrZSBhbnkgaW50ZXJ2YWwgd2l0aCBhIG1lYW4gYWJvdmUgYWJvdXQgMi4wIGtXaCBzaG91bGQgcXVhbGlmeSBhcyBhbiBFViBjaGFyZ2luZyBzZXNzaW9uLiAKLSBEbyB3ZSBuZWVkIHRvIGRvIGFueXRoaW5nIHRvIGF2b2lkIG1peHVwcyB3aXRoIGhpZ2ggSFZBQyB1c2FnZT8gCi0gRG8gd2UgbmVlZCB0byBwdXQgYW4gaWRlYWwgbGVuZ3RoIG9mIGNoYXJnZSBvbiB0aGUgaW50ZXJ2YWxzPyBNYXhpbXVtIHZhcmlhdGlvbiBkdXJpbmcgdGhlIHBlcmlvZD8gCi0gVXNlIHBvc3QucHJvYiB0byBmaWx0ZXIgb3V0IGFueSBpbnRlcnZhbHM/IAoKYGBge3J9CnJlc3VsdCA8LSBzdW1tYXJ5KGJjcCwgbm9ybS50ZXN0ID0gVFJVRSkKCmJjcF9zZWdtZW50cyA8LSBhc190aWJibGUocmVzdWx0JHJlc3VsdCRzZWdtZW50KSAlPiUgCiAgbXV0YXRlKGxlbmd0aCA9IGVuZCAtIGJlZ2luKSAlPiUgCiAgbXV0YXRlKGxhZ19tZWFuID0gbGFnKG1lYW4sMSkpICU+JSAKICBtdXRhdGUobGVhZF9tZWFuID0gbGVhZChtZWFuLDEpKSAlPiUgCiAgbXV0YXRlKGFkaiA9IChsYWdfbWVhbitsZWFkX21lYW4pLzIpCgpiY3Bfc2VnbWVudHMgJT4lIAogIGZpbHRlcihtZWFuID4gMi4wICYgbWVhbiA8IDMuMCkgCgpgYGAKCmBgYHtyfQoKc2FtcGxlX2FkaiA8LSAgc2FtcGxlICU+JSAgCiAgbGVmdF9qb2luKGJjcF9zZWdtZW50cywgYnkgPSBqb2luX2J5KHRpbWVfeD49YmVnaW4sdGltZV94PD1lbmQpKSAlPiUgCiAgbXV0YXRlKGNvbnN1bWVkX2tXaF9hZGogPSBjYXNlX3doZW4oCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbWVhbiA+IDIuMCAmIGxlbmd0aCA+PSA0ICYgU0QgPCAxLjEgfiBjb25zdW1lZF9rV2ggLSAobWVhbiAtIGFkaiksIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIC5kZWZhdWx0ID0gY29uc3VtZWRfa1doCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgKQogICAgICAgICApICU+JSAKICBtdXRhdGUoY29uc3VtZWRfa1doX2FkaiA9IGNhc2Vfd2hlbigKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb25zdW1lZF9rV2hfYWRqIDwgMCB+IDAsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIC5kZWZhdWx0ID0gY29uc3VtZWRfa1doX2FkagogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICkKICApCgpzYW1wbGVfYWRqICU+JSAKICBtdXRhdGUoZGlmZiA9IGNvbnN1bWVkX2tXaCAtIGNvbnN1bWVkX2tXaF9hZGopICU+JSAKICBncm91cF9ieShzb2xhcl95ZWFyKSAlPiUgCiAgc3VtbWFyaXplKGRpZmYgPSBzdW0oZGlmZikpCgpnZ3Bsb3Qoc2FtcGxlX2FkaiwgYWVzKGRhdGV0aW1lLCBjb25zdW1lZF9rV2hfYWRqLCBncm91cD1zb2xhcl95ZWFyLCBjb2xvcj1zb2xhcl95ZWFyKSkgKwogICNmYWNldF9ncmlkKHJvd3MgPSB2YXJzKG1vbnRoKSkgKwogIGdlb21fcG9pbnQoc2l6ZT0wLjUpICsKICBnZ3RpdGxlKCJIb21lIEVsZWN0cmljaXR5IGluIDE1LU1pbnV0ZSBJbnRlcnZhbHMiKSArCiAgeWxhYigiRWxlY3RyaWNpdHkgQ29uc3VtcHRpb24gKGtXaCkiKQoKYGBgCgoKc29sYXJfeWVhcgo8ZmN0cj4KZGlmZgo8ZGJsPgowCS0wLjM2NzAJCQkKMQk1OTUuODc3MAkJCQoyCTM1MDQuMTM4NwkJCQozCTMyMS45NDIxCQkJCgp2cy4gCgpzb2xhcl95ZWFyCjxmY3RyPgprV2gKPGRibD4KMQk2ODMJCQkKMgk0MzgyCQkJCjMJNjA3CgpVbmRlcmVzdGltYXRlcyAgCgpJZ25vcmUgdGhlIGJjcCBtb2RlbCwganVzdCBwdWxsIGRhdGEgb24gaG9tZSBjaGFyZ2luZyBmcm9tIEVtcG9yaWEKCmBgYHtyfQplbXBvcmlhX2RhaWx5IDwtIHJlYWRfY3N2KCcvVXNlcnMvYmVuamVkbG92ZWMvSW9uaXFDaGFyZ2luZy9lbXBvcmlhX2RhaWx5LmNzdicpICU+JSAKICBtdXRhdGUoZGF0ZSA9IGRhdGUoZGF0ZXRpbWUpKQoKZGFpbHlfZWxlY3RyaWNpdHlfYWRqIDwtIGRhaWx5X2VsZWN0cmljaXR5ICU+JSAKICBtdXRhdGUoZGF0ZSA9IGFzX2RhdGUoZGF0ZSkpICU+JSAKICBsZWZ0X2pvaW4oZW1wb3JpYV9kYWlseSwgYnkgPSBqb2luX2J5KGRhdGUpKSAlPiUgCiAgbXV0YXRlKGNvbnN1bWVkX2tXaF9hZGogPSBjb25zdW1lZF9rV2ggLSByZXBsYWNlX25hKGtXaCwwKSkgJT4lIAogIHJlbmFtZShzb2xhcl95ZWFyID0gc29sYXJfeWVhci54LCBjaGFyZ2luZ19rV2ggPSBrV2gpCgpkYWlseV9lbGVjdHJpY2l0eV9hZGogJT4lIAogIG11dGF0ZShkaWZmID0gY29uc3VtZWRfa1doIC0gY29uc3VtZWRfa1doX2FkaikgJT4lIAogIGdyb3VwX2J5KHNvbGFyX3llYXIpICU+JSAKICBzdW1tYXJpemUoY29uc3VtZWRfa1doID0gc3VtKGNvbnN1bWVkX2tXaCksIGNoYXJnaW5nX2tXaCA9IHN1bShyZXBsYWNlX25hKGNoYXJnaW5nX2tXaCwwKSksIGNvbnN1bWVkX2tXaF9hZGogPSBzdW0ocmVwbGFjZV9uYShjb25zdW1lZF9rV2hfYWRqLDApKSwgZGlmZiA9IHN1bShyZXBsYWNlX25hKGRpZmYsMCkpKQoKZ2dwbG90KGRhaWx5X2VsZWN0cmljaXR5X2FkaiwgYWVzKHNvbGFyX2RheSwgY29uc3VtZWRfa1doX2FkaiwgZ3JvdXAgPSBzb2xhcl95ZWFyLCBjb2xvciA9IHNvbGFyX3llYXIpKSArIAogIGdlb21fcG9pbnQoKSArCiAgI3NjYWxlX3hfY29udGludW91cyhicmVha3MgPSBjKDEsIDMyLCA2MCwgOTEsIDEyMSwgMTUyLCAxODIsIDIxMywgMjQ0LCAyNzQsIDMwNSwgMzM1KSwgbGFiZWxzID0gbW9udGguYWJiKSArCiAgdGhlbWUoYXhpcy50ZXh0LnggPSBlbGVtZW50X3RleHQoYW5nbGUgPSA5MCkpIAoKCmBgYAoKYGBge3J9Cm1vbnRobHlfcmVkdWN0aW9uIDwtIGRhaWx5X2VsZWN0cmljaXR5X2FkaiAlPiUgCiAgbXV0YXRlKHllYXJfMSA9IHN1bShjYXNlX3doZW4oc29sYXJfeWVhciA9PSAxIH4gY29uc3VtZWRfa1doX2FkaikpLCB5ZWFyXzIgPSBzdW0oY2FzZV93aGVuKHNvbGFyX3llYXIgPT0gMiB+IGNvbnN1bWVkX2tXaF9hZGopKSkgJT4lIAogIGdyb3VwX2J5KG1vbnRoID0gbW9udGgoZGF0ZSkpICU+JSAKICBzdW1tYXJpemUoeWVhcl8xID0gcm91bmQoc3VtKHJlcGxhY2VfbmEoeWVhcl8xLDApKSksIHllYXJfMiA9IHJvdW5kKHN1bShyZXBsYWNlX25hKHllYXJfMiwwKSkpKSAlPiUgCiAgbXV0YXRlKHBjdCA9IHJvdW5kKCh5ZWFyXzEteWVhcl8yKS95ZWFyXzEsMikpCgptb250aGx5X3JlZHVjdGlvbgoKZGFpbHlfZWxlY3RyaWNpdHlfYWRqICU+JSAKICBncm91cF9ieShtb250aCA9IG1vbnRoKGRhdGUpLCBzb2xhcl95ZWFyKSAlPiUgCiAgc3VtbWFyaXplKGNvbnN1bWVkX2tXaF9hZGogPSBzdW0oY29uc3VtZWRfa1doX2FkaikpCgpgYGAKCg==